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By means of the first-principles calculations combined with the tight-binding approximation, the 
strain-induced semiconductor-semimetal transition in graphdiyne is discovered. It is shown that 
the band gap of graphdiyne increases from 0.47 eV to 1.39 eV with increasing the biaxial tensile 
strain, while the band gap decreases from 0.47 eV to nearly zero with increasing the uniaxial tensile 
strain, and Dirac cone-like electronic structures are observed. The uniaxial strain-induced changes 
of the electronic structures of graphdiyne come from the breaking of geometrical symmetry that lifts 
the degeneracy of energy bands. The properties of graphdiyne under strains are disclosed different 
remarkably from that of graphene. 
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Carbon, one of the most common chemical elements 
in Nature, has a lot of allotropes, plenteous properties 
and numerous potential applications. Besides the well- 
known graphite, graphene 1 , and diamond, many other 
novel carbon allotropes, such as T-carbon 2 -, M-carbon 3 , 
graphyne 4 ^, graphdiyne 6 , octgraphene 7 , etc., are pro- 
posed and studied. Owing to its unique electronic 8 , 
mechanical and thermal properties, graphene stimulated 
considerable interest in itself and other two-dimensional 
(2D) materials 9 12 . Graphene's brothers, graphyne and 
graphdiyne, are also intriguing 2D carbon materials. In 
particular, both building blocks and cut-outs were al- 
ready obtained experimentally^ 3 - - — . It is known that 
graphene is a zero-gap semimetal, while graphdiyne is 
a semiconductor with a direct band gap. In addition, 
graphdiyne has the high-temperature stability and shows 
mechanical properties similar to graphene^, which is 
considered to be a potential material for nanoelectron- 
ics. The electronic modulation of graphdiyne^"— is thus 
of great interest. 

As graphene has no energy gap, which poses diffi- 
culties for wide applications in microelectronic devices, 
plenty of researches were devoted to open a band gap 
in graphene. One method, among others, is to invoke 
strains 22 28 , which is also found effective in modulat- 
ing electronic properties of graphdiyne and other 2D 
materials such as BN shee t 29 ! 30 . In this work, a sys- 
tematic first-principles density-functional theory (DFT) 
study on the strain-induced changes of electronic struc- 
tures in graphdiyne was performed. In the absence of 
a strain graphdiyne is known as a semiconductor with 
a direct band gap. In the presence of strains, we found 
that the band gap of graphdiyne increases with increas- 
ing the biaxial tensile strains, but the band gap decreases 
with increasing a uniaxial tensile strain either along the 
armchair or zigzag direction. The electronic properties 
around the Fermi surface are dominated by 2p z orbit als, 
and the Dirac cone-like electronic structures can be ob- 




FIG. 1. (Color online) The schematic structures of graphdiyne 
under various tensile strains, (a) under biaxial tensile strains, 
where the parallelogram (dot line) indicates a unit cell with 
length L x and width L y ; (b) under a uniaxial strain along 
the armchair direction; (c) under a uniaxial strain along the 
zigzag direction. The corresponding Brillouin zones are illus- 
trated below each structure. 



tained. The calculated electron density of graphdiyne un- 
der different strains shows that C-C bonds are elongated, 
on which the electron densities are reduced significantly, 
while C = C bonds are changed little. We also proposed 
a tight-binding model of tt electrons to describe the low- 
energy physics around the Fermi level, which gives the 
results in well agreement with the DFT calculations. The 
decline of the band gap of graphdiyne under an asymmet- 
ric, uniaxial strain may be owing to the breaking of the 
geometrical symmetry that lifts the degeneracy of energy 
bands. These observations in graphdiyne are different 
from those of graphene, and may be useful for nanoelec- 
tronics. 

All calculations were performed by DF T 31 i 32 imple- 
mented within the Vienna ab initio simulation pack- 
age (VASP) 33,34 with the projector augmented-wave 
(PAW) method 3 -^, and the results are rechecked using the 
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FIG. 2. (Color online) (a) L y versus L x and (b) the band 
gap versus L x of graphdiyne under biaxial tensile strains; (c) 
L y versus L x and (d) the band gap versus L x of graphdiyne 
under a uniaxial strain along the armchair direction; (e) L x 
versus L y and (f) the band gap versus L y of graphdiyne under 
a uniaxial strain along the zigzag direction. The strain- stress 
relation for each case is shown in the inset of (a), (c) and (e), 
respectively. The crosses (red) indicate the largest strain we 
considered for each case, which can be taken as the transition 
point from elastic to nonelastic situation. 



Quantum-ESPRESSO package^. The generalized gra- 
dient approximation (GGA) developed by Perdew and 
Wang 37 are adopted for the exchange correlation poten- 
tial. The Monkhorst-Pack scheme is used to sample the 
Brillouin zone^, and a mesh ofl7x 17 x 1 /c-point sam- 
pling was used for the calculations. The total energy was 
converged to within 1 meV with the plane-wave cutoff 
energy 500 eV. The geometries were optimized when the 
remanent Hellmann-Feynman forces on the ions are less 
than 0.01 eV/nm. The distances between two layers in 
the super cell are more than 10 A to avoid the inter layer 
interactions. 

The schematic structure of graphdiyne is shown in Fig. 
1(a), which can be considered as the structure that one 
third of C-C bonds in graphene are displaced with two 
C=C units. These two C=C units form acetylenic link- 
ages. The parallelogram in Fig. 1(a) indicates the unit 
cell of graphdiyne with length L x and width L y . The ge- 
ometrical and electronic properties of graphdiyne in the 
absence of strains were calculated, and the results re- 
veal that the lattice constants are L x = 0.945 nm and 
Ly = 0.82 nm, and the a direct band gap is 0.47 eV at Y 
point, as shown in Fig. 3(a), which are in good agreement 
with the previous calculations—^. 



Now let us consider graphdiyne under three types of 
strains: (i) under symmetrical biaxial tensile strains, as 
indicated in Fig. 1(a); (ii) under a uniaxial tensile strain 
along the armchair direction, as indicated in Fig. 1(b); 
and (iii) under a uniaxial tensile strain along the zigzag 
direction, as indicated in Fig. 1(c). For the case (i), when 
graphdiyne is stretched under biaxial strains, L x and L y 
are stretched uniformly in the same ratio, and the strain 
can be represented as a ratio defined by (L x — L x )/L x 
in percent. For the cases with uniaxial strains, only one 
direction is stretched, while another direction is set to 
free, which makes the whole system reach the lowest en- 
ergy. Thus, for case (ii) the strain can be measured by 
(L x — L x )/L x in percent, and for case (iii), the strain 
can be measured by (L y — Ly)/Ly in percent. In our 
calculations, the strain is changed in every 1%. 

Fig. 2(a) shows L y versus L x and stress-strain relation 
of graphdiyne under symmetrical biaxial strains. L y vs. 
L x curve gives a raising straight line, and the stress-strain 
curve shows a minimum at the strain of 15%, which can 
be viewed as the transition point from elastic to nonelas- 
tic situation for graphdiyne under the biaxial strains. In 
this case, the band gap increases from 0.47 eV to 1.39 
eV with increasing L x , as seen in Fig. 2(b). When the 
strain is executed along the armchair direction (case (ii)), 
L y vs. L x shows a declining straight line, as revealed in 
Fig. 2(c), which is similar to that of graphene under a 
strain perpendicular to C-C bonds (Ref. [23), indicating 
that graphdiyne is the 2D material with a positive Pois- 
son's ratio. The elastic-nonelastic transition occurs at 
the strain of 16%, close to the value of graphdiyne under 
biaxial strains. In this case, the band gap of graphdiyne 
is decreased almost linearly from 0.47 to 0.02 eV with 
increasing L x (Fig. 2(d)), which is distinct from the case 
under biaxial strains. For case (iii), L y as a function of 
L x shows a concave descending curve, as shown in Fig. 
2(e), which is also similar to that of graphene under a 
strain parallel to C-C bonds (Ref. l23l ). The elastic- 
nonelastic transition occurs around the strain of 20%. 
Fig. 2(f) presents the band gap of graphdiyne under the 
strain along the zigzag direction, which decreases from 
0.47 to 0.02 eV with increasing L y from 0.82 nm to 0.88 
nm. When L y is between 0.90 nm to 0.95 nm, the band 
gap is nearly zero, indicating a semiconductor-semimetal 
transition. When L y is larger than 0.95 nm, the band gap 
re-opens again. This latter property is also supported by 
the TB calculations (see below). These properties are 
dramatically different from those of graphen o 23 i 41 i 42 . 

In order to better understand the strain effect on elec- 
tronic properties of graphdiyne, the energy bands, partial 
density of states (PDOS), and electron density distribu- 
tion under biaxial or uniaxial strains are calculated. Fig. 
3(a) shows the electronic structure of graphdiyne in the 
absence of strain. It is seen that graphdiyne is a semicon- 
ductor with a direct band gap 0.47 eV at Y point, which is 
consistent with previous result a 39 i 4Q . The degeneracy of 
energy bands can be observed in both the lowest conduc- 
tion bands and the highest valence bands near Y point. 
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FIG. 3. (Color online) The band structures calculated by the 
density functional theory (DFT) and the tight-binding ap- 
proximation (TBA), the partial density of states (PDOS) and 
the three-dimensional energy bands at Y point for graphdiyne 
(a) without a strain, and under a symmetrical biaxial tensile 
strain of (b) 5%, (c) 9%, and (d) 15%, respectively. 



The PDOS calculations show that 2p z orbitals dominate 
the states around the Fermi level, while 2p x and 2p y or- 
bitals contribute the same, which are mainly located at 
least 1 eV below the Fermi level. The three-dimensional 
(3D) energy bands near Y point around the Fermi level 
reveal two opposite bowls, showing the parabolic-like dis- 
persions. Figs. 3(b)-(d) show the electronic structure of 
graphdiyne under a symmetrical biaxial tensile strain of 
5%, 9% and 15%, respectively. It can be observed that 
for all these cases, the minimal energy gap is still lo- 
cated at T point, but it gradually becomes larger and 
larger as the strain increases. The degeneracy of the low- 
est conduction bands and the highest valence bands still 
maintains, as the geometrical symmetry is unchanged un- 
der the symmetrical strain. Despite of the variation of 
band gaps, the PDOS remains almost unchanged, and 
2p z orbitals are still dominating around the Fermi level. 
From these figures, one may note that the energy bands 



> 

CD 



(b) 



1 


f ^ 


J2s. 

" — TBA 


' — -p, 


'jsL 















-0.4 

lfl.2 

| -0.2 
'-Q.4 



(c) 



TBA _. f 



-0.4 
0.2 


-0.2 
'■0.4 



> 
CD 




FIG. 4. (Color online) The band structures calculated by the 
density functional theory (DFT) and the tight-binding ap- 
proximation (TBA), the partial density of states (PDOS) and 
the three-dimensional energy bands at Y point for graphdiyne 
(a) without a strain, and under a uniaxial tensile strain along 
the armchair direction with the strain of (b) 5%, (c) 9% and 
(d) 15%, respectively. 



become more flat as the strain increases, revealing that 
7r electrons become more localized under the symmetri- 
cal biaxial strains. This observation is also confirmed in 
the electron densities shown in Figs. 6(a)-(c), where the 
electron densities between carbon atoms decrease signif- 
icantly. 

Figures 4(a)- (c) present the electronic structures of 
graphdiyne under uniaxial tensile strains along the arm- 
chair direction with the strain of 5%, 9% and 15%, re- 
spectively. It is found that the minimal energy gaps are 
still located at Y point. However, with the increase of 
strain, the degeneracies of the highest valence bands and 
the lowest conduction bands are lifted. The energy levels 
of the highest valence bands move up and those of the 
lowest conduction bands move down, resulting in that the 
energy gaps diminish with the increase of strain. Espe- 
cially, the energy gap becomes very small (0.02 eV) when 
the strain is larger than 15%. A careful observation finds 
that the highest valence bands and the lowest conduction 
bands from Y point to M point are separated much more 
than the bands from Y point to S point, and the gap at M 
point decreases gradually, which is resulted from the bro- 
ken symmetry in the Brillouin zone. It is also uncovered 
in PDOS calculations that 2p z orbitals still dominate the 
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FIG. 5. (Color online) The band structures calculated by 
the density functional theory (DFT) and the tight-binding 
approximation (TBA), the partial density of states (PDOS) 
and the three-dimensional energy bands for graphdiyne under 
a uniaxial tensile strain along the zigzag direction with the 
strain of (a) 5% and (b) 9% at F point; and (c) 15% and (d) 
19% at the k point with minimal energy gap. 

states around the Fermi level, and the contribution of 
2p x and 2p y orbitals becomes different. In addition, it 
is interesting to observe that a Dirac cone-like electronic 
structure is formed (Fig. 4(c)). 

The electronic structures of graphdiyne under uniaxial 
tensile strains along the zigzag direction are more com- 
plex, as displayed in Figs. 5(a)-(d), in which the cases 
with strains of 5%, 9%, 15% and 19% are presented. 
When the strain is increased from zero to 9% (Figs. 5(a) 
and 5(b)), the degeneracy of the highest valence bands 
is removed, and the minimal band gap opens at T point 
and tends to decrease to nearly zero (Fig. 5(b)). When 
the uniaxial strains between 9% and 17% are applied, 
the band gap still maintains nearly zero but the posi- 
tion at which the gap closes is moving gradually from T 
point to S point. The corresponding 3D band structure 
around the Fermi level manifests a Dirac cone-like struc- 




FIG. 6. (Color online) The distributions of electron densities 
of graphdiyne under (a) no strain; biaxial strains of (b) 9% 
and (c) 15%; uniaxial tensile strains along the armchair di- 
rection with the strain of (d) 5%, (e) 9% and (f) 15%, and 
along the zigzag direction with the strain of (g) 5%, (h) 15% 
and (i) 19%. 



ture, in which the cone becomes sharp with the increase 
of strains, as seen in Fig. 5(c). This observation shows 
that a strain-induced semiconductor-semimetal transi- 
tion takes place in graphdiyne. When the strain is larger 
than 17%, the band gap opens again at S point (Fig. 
5(d)). It is interesting to note that, in contrast to the 
case for the strain along the armchair direction, in this 
case the bands near the Fermi level from T to S point 
are lifted much more than the bands from T to M point, 
which is due to the broken symmetry of the Brillouin 
zone. The PDOS calculations reveal that 2p z orbitals 
dominate the states around Fermi level, while 2p x and 
2p y orbitals contribute in different ways, which is also 
owing to the symmetry breaking caused by the asym- 
metrical uniaxial tensile strains. 

It is known that graphdiyne is composed of acetylenic 
linkages (two linked C=C bonds) and carbon hexagons. 
The distributions of electron densities of graphdiyne un- 
der different tensile strains are given in Fig. 6. The 
distribution of electron density of graphdiyne without a 
strain is shown in Fig. 6(a). As one may see, the elec- 
tron densities on carbon hexagons are similar to those 
of graphene, which display a six- fold symmetry. For 
acetylenic linkages, the electron density of C=C bonds is 
obviously higher than that of the C-C bonds. As shown 
in Figs. 6(b) and (c), with increasing the symmetrical 
biaxial tensile strain, the electron densities on carbon 
hexagons keep the symmetry, but the densities between 
carbon atoms become smaller. In acetylenic linkages, 
the electron densities on C=C bonds keep unchanged, 
while those of C-C bonds decrease significantly, showing 
that C-C bonds become weak on carbon hexagons and 
acetylenic linkages, which may cause the energy bands 
more narrow and fiat, leading to the band gaps become 
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TABLE I. The lattice constant, plane density p, energy gap E g between the bottom of conduction band and the top of valence 
band, Poission's ratio v, and the angle shown in Fig. 1(a) under different strains for graphdiyne. 



strain% 


lattice constant (A) 


p (mg/m 2 ) 


E g (eV) 


V 


6 (°) 


E c (eV/a 


unstrain (present work) 


9.45 


0.466 


0.47 


0.38, 0.476 


120.00 


8.48 


unstrain (Ref. ^ and 40) 


9.44 


0.466 


0.53 


0.41, 0.44 


120.00 


7.78 


biaxial (5%) 


9.823 


0.426 


0.677 




120.00 


8.34 


biaxial (9%) 


10.30 


0.396 


0.888 




120.00 


8.10 


biaxial (15%) 


10.87 


0.352 


1.389 




120.00 


7.58 


armchair direction (5%) 


9.397 


0.458 


0.321 


0.386 


111.616 


8.45 


armchair direction (9%) 


9.374 


0.449 


0.156 


0.384 


108.213 


8.37 


armchair direction (15%) 


9.366 


0.437 


0.037 


0.377 


107.255 


8.20 


zigzag direction (5%) 


9.755 


0.456 


0.294 


0.434 


127.397 


8.46 


zigzag direction (9%) 


10.022 


0.445 


0.024 


0.387 


131.822 


8.36 


zigzag direction (15%) 


10.435 


0.427 


0.019 


0.335 


136.836 


8.15 


zigzag direction (19%) 


10.717 


0.421 


0.299 


0.308 


140.198 


7.99 



larger. The electron densities of graphdiyne under the 
uniaxial tensile strain along the armchair direction are 
presented in Figs. 6(d)-(f). As the carbon linkages along 
the armchair direction are elongated, the electron densi- 
ties on C-C bonds decrease while those of C=C bonds 
change little. The carbon linkages along other directions 
are compressed, where the electron densities on C-C and 
C=C bonds are still high and even slightly larger than 
those of the unstrained situation. Figs. 6(g) and (h) re- 
veal the electron densities of graphdiyne under the uni- 
axial tensile strain of 5% and 15% along the zigzag di- 
rection, respectively. In contrast to the case along the 
armchair direction, the carbon linkages along the arm- 
chair direction are compressed, and the electron densities 
are slightly larger than those of the unstrained situation, 
while the carbon linkages along the zigzag direction are 
elongated, leading to the electron densities on C-C bonds 
decrease. However, when the strain along the zigzag di- 
rection is applied up to 19% (Fig. 6(i)), the electron 
densities on carbon hexagons and C-C bonds decrease 
dramatically, and a large number of electrons are local- 
ized on C=C bonds, which may be the reason that the 
energy gap of graphdiyne opens again. 

To understand the DFT results of the electronic prop- 
erties of graphdiyne under different tensile strains, a 
tight-binding (TB) model of tt electrons in graphdiyne 
is proposed near the Fermi level. As indicated above in 
PDOS analyses, the electronic states around the Fermi 
level are predominantly contributed by 2p z orbitals of 
carbons in graphdiyne, no matter whether strains are ap- 
plied or not. Thus, only tt electrons and nearest neighbor 
hoppings are sufficient to describe the low-lying proper- 
ties of electrons in graphdiyne. The TB Hamiltonian 
could be written as 

H=Y / tiA& j + h.c.) (1) 

(iJ) 

where the operators c\ and q are the creation and anni- 




FIG. 7. (Color online) The atom sites i (from 1 to 18) in 
one unit cell (the parallelogram) of graphdiyne are labeled 
with numbers, and Uj in Eq. (1) is defined as the hopping 
amplitude between sites i and j. 



hilation operators of tt electrons at site z, denotes 
the sum over nearest neighbors within the unit cell, and 
stands for the hopping amplitude between sites i and 

j- 

As shown in Fig. 7, there are 18 atoms in one unit 
cell. As i and j should go through all of these atoms, the 
matrix form of Hamiltonian is 18x18, where the number 
of independent hopping amplitude Uj relies on the sym- 
metry of the unit cell. We first calculate the TB band 
structure of graphdiyne without a strain, as given in Fig. 
3(a). One may see that the TB results are in good agree- 
ment with the DFT calculations, implying that the low- 
energy properties of tt electrons in grpahdiyne can be well 
described by the TB Hamiltonian with four independent 
hopping amplitude tij. When the biaxial tensile strain is 
applied to graphdiyne, the corresponding TB band struc- 
tures are found very well consistent with the DFT results 
(Figs. 3(b)-3(d)), where there is also 4 independent tij 
parameters owing to the same symmetry as in the primi- 
tive graphdiyne. We also found that the fitting values of 
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tij become smaller with increasing the strain, and tij on 
C-C bonds are much smaller than those on C=C bonds, 
which are also consistent with the decrease of electron 
densities shown in Figs. 6(b) and (c). 

For the cases of graphdiyne under uniaxial tensile 
strains, there are eight independent parameters Uj. The 
corresponding TB band structures are shown in Figs. 4 
and 5, which give the results in well agreement with the 
DFT calculations. In these cases the degeneracies of the 
lowest conduction bands and the highest valence bands 
are lifted, and the band gaps become narrow. For the 
case of strains applied along the armchair direction, Uj 
along this direction are smaller than those along the other 
directions, being consistent with the electron densities 
shown in Figs. 6(d)-(f). The minimal band gap is con- 
firmed at T point, and decreases to nearly zero. The case 
of tensile strains applied along the zigzag direction also 
recovers the DFT results (Figs. 6(g)-(i)). The TB calcu- 
lations support as well that the broken geometrical sym- 
metry of graphdiyne caused by the asymmetrical uniaxial 
strains lifts the degeneracy of energy bands, resulting in 
the shrink of the band gap. 

By means of first-principles calculations combined with 



the TB approximations, the effect of symmetrical biaxial 
and asymmetrical uniaxial strains on electronic proper- 
ties of graphdiyne is explored, which is shown remark- 
ably different from that on graphene. The calculated 
results reveal that the symmetrical biaxial strain can en- 
large the band gap of graphdiyne, while the asymmetri- 
cal uniaxial strains could remove the band degeneracies, 
giving rise to the decrease of the band gap. It is found 
that the electronic structures around the Fermi level are 
dominantly determined by 2p z orbit als of tt electrons. 
The strain-induced semiconductor-semimetal transition 
in graphdiyne is observed when the asymmetrical uniax- 
ial tensile strain is applied, and Dirac cone-like electronic 
structures are obtained. This study provides a possible 
way to modulate the electronic properties of 2D materi- 
als like graphdiyne through external strains, which may 
have potential applications. 
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